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Abstract 

Although various physical systems have been explored to produce entangled states 
involving electromagnetic, strong and weak interactions, the gravity has not yet been 
touched in practical entanglement generation. Here, we propose an experimentally 
feasible scheme with current technology for generating spin entangled neutron pairs via 
the Earth's gravitational field. The scheme is realized by passing two neutrons through 
a specific rectangular cavity, where the gravity adjusts the neutrons into entangled 
state. This provides a simple and practical way for the implementation of the test of 
quantum nonlocality and statistics in gravitational field. 

Ever since the discovery of Bell inequalities PQ, the generation of entanglement with 
various physical systems has been the intensively studied subject. Now the entanglement 
can be generated not only in optical [2j [3] , atomic [4], solid state [5] systems where only the 
electromagnetic interactions is involved, but also with baryons [6], leptons [7J, and mesons 
[3 El HOI HI] where the strong or weak interaction emerges as the dominant force. Many 
of these entanglement generation schemes have been experimentally realized in testing the 
violation of Bell's inequalities which reveal the nonlocality of quantum theories. Some of 
them have further found their roles in quantum computations and quantum information 
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processing [121 021 EH] • One considerable interest of exploring these various entangled systems 
is to show that the nonlocal correlation is not a peculiarity attributed to specific interactions 
but a universal quantum phenomena. Besides, such practical entangled systems contribute 
also as the crucial physical resource for quantum information science. 

It has been noticed that three of the four fundamental interactions in nature, i.e., elec- 
tromagnetic, strong and weak, are capable of generating entanglement leaving the gravity 
a sole exception [9J- Although some quantum effects of the classical gravity have been ob- 
served, i.e., quantum interference [15], discrete energy levels of neutrons in the gravitational 
potential [16], there is still no report on quantum entanglement generation concerning the 
Earth's gravitational filed. In the above two cases, the neutrons act as a messenger revealing 
quantum effects of gravity. Moreover, the neutrons have long been considered the ideal tool 
to test the Bell inequalities [17] . However, due to the extremely small neutron-neutron scat- 
tering length (~ 10~ 14 m [H]), the usual way to generate entangled neutron pairs through 
their low energy scattering is a considerable technical challenge. Up to now, only the en- 
tanglement among different degree of freedoms of a single neutron has been experimentally 
realized [HI [20]. 

Here, we propose a scheme to generate the Einstein-Podolsky-Rosen state via Earth's 
gravitational field. The scheme composed of three different functional components: an 
energy filter that monochromatizes the neutrons, a rectangular cavity which entangles two 
neutrons, and the neutron spin polarization analyzers. The main idea is to guide a pair of 
monochromatic neutrons into the cavity where they are enforced into the same energy state. 
The spin singlet state should then be obtained due to the antisymmetrization requirement 
of the identical fermions. We design a practical experimental setup for our scheme, by which 
we show that, with the current technology, the gravity dominant entanglement generation 
and the nonlocality test with entangled neutron pairs are within the experimental reach. 

As has been observed in [IB], neutrons falling towards a horizontal reflecting mirror will 
distributed discontinuously in the vertical direction. Such a system can be described by the 
quantum theory of a particle bouncing in the gravitational field above a perfect mirror |21j. 
The Schrodinger equation governing the motion of the neutrons in the vertical dimension 
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Figure 1: A rectangular cavity where the two lower surfaces are neutron mirrors and the two 
upper surfaces are neutron absorbers, (a): A neutron N passes through the rectangular cavity 
with longitudinal velocity v z ; (b): The transection of the rectangular cavity, h is the height of the 
neutrons relative to the zero gravitational potential surface. 



reads 



2m d£ 2 




Here £ is the height of neutron from the horizontal mirror. Equation ([T]) can be solved and 
shows discrete energy eigenstates 



(2m 2 a) 1/3 2 1 / 3 
M€) = ^nAi(eL_|L_ _ £ n ____) , (2 ) 

(mg 2 ^ 2 ) 1 / 3 

where Ai are Airy functions and Af n is a normalization constant, E n = a n — jyjj cor ~ 

responds to the energy eigenvalue of state <p n with a n being the nth zero of Airy function 
Ai(-cv) = 0. 

Now considering a specifically designed rectangular cavity of Fig. [U where the two lower 
surfaces of the cavity are neutron mirrors while the two upper surfaces are neutron absorbers. 
For each neutron in this cavity, it is subjected to a constant gravitational force F = mg, 
where m is the mass of neutron, g is the acceleration constant near the Earth's surface. To- 
gether with the lower mirrors, the Earth's gravitational field provides the confining potential 
well for the neutrons, which is 

mg(x + y)/ \/2 x > and y > 

VU.y)- { ' , (3) 

-foo x < or y < 
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where we have chosen the zero potential surface in Fig. [TJ In the transection of the cavity, the 
X-Y plane of Fig. Hfb), the Schrodinger equation of motion for neutron takes the following 
form 

h 2 ( d 2 d 2 \ 

+ —)^(x,y) + V(x,y)^(x,y)=E^(x,y). (4) 



2m \dx 2 dy 2 / 

The energy eigenstates and eigenvalues of this equation can be similarly obtained as that of 
equation (JTJ), and can be simply formulated as 

^n, m (x, y) = ip n (x)ip m (y) , E n>m = E n + E m . (5) 

Here, ip n (x) = J\f n Ai(x/l — E n /e ), Af n is the normalization constant; l and e are the 
characteristic length and energy defined as 

l = h 2/3 /( y V2m 2 g) 1/3 ~ 6.59 • KT 6 m , (6) 



e = {/mg 2 h 2 /4: ~ 4.78 • 1(T 13 eV ; (7) 

and E n = a n e is the eigen energy of tj) n with a n being defined after equation fl2]). The first 
four lowest eigenstates are ^0,0) ^o,i ; ^i,cb an d ^1,1, with \^ n>m {x, y)\ 2 being the value of the 
neutron distribution densities in the X-Y plane, see Fig. |2j Their corresponding energies 
can be listed as follows 

£ ,o = E + E = a e + a e ~ 2.22 x 10" 12 eV , (8) 

Eo tl = Eo + E l = a 6 + aie ~ 3.06 x 10~ 12 eV , (9) 

E 1>0 = E 1 +E = ai£ + a e ~ 3.06 x 10" 12 eV , (10) 

E ltl = E x + E 1 = a^o = ai£ ~ 3.90 x 10" 12 eV , (11) 

where Eo,o is the ground state energy and E 0t i and E^ are energies of the first degenerated 
excited states. 

From the density distributions \^f n ,m(x, y)\ 2 plotted in Fig. [2j we see that the quantum 
states of higher energies become densely distributed in the area with larger values of x, y. If 
the two upper absorbers are set at a position of x — y ~ 31q, the transection of rectangular 
cavity would have the size holding only the ground state. In this configuration, where the 
excited states ^0,1, ^1,0, ^1,1 and other even higher excited states are substantially absorbed 





Figure 2: The probability distributions of the first four lowest eigenstates in the transection of 
the rectangular cavity: |^o,o| 2 ; l^o.il 2 ; |^i,o| 2 j |^i,i| 2 - The X and Y axes are in unit of Iq. 

while the neutron transmitting in the cavity, only the ground state may survive. The minimal 
time At while the neutron has to stay in the cavity to resolve different quantum states is 
[22] 

At = h/(E 0)1 - E 0)0 ) ~ 7.84 • l(r 4 s . (12) 

After a significantly longer time duration T ^ At in the cavity, a substantial absorbtion of 
the excited state will be expected in the present configuration of x = y ~ 3/o- 

The neutrons from the same monochromatic source having the the same energy spectrum 
and velocity will have the same uncertainty in momentum, the Ap z . The Heisenberg uncer- 
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tainty relation ascertains the indistinguishability of the neutron in longitudinal position z 
as: AzAp z > h/2. For slow neutrons Ap z = mAv z , we can define the full width at half max- 
imum (FWHM) of the velocity spectrum Av z whose relation with Av z can be estimated via 
the minimum uncertainty wave packet: Av z = 2\/2 m~2- Av z . Thus we have Az Av z > H/m. 
Two neutrons A and B will keep a relative displacement of 5z = T5v z = T\vz — vi \ 
through the passage in the cavity. In case the relative displacement of neutrons 5z < Az, 
i.e., Sz < h/(mAv z ), while they enter the cavity "simultaneously", two neutrons will stay 
longitudinally indistinguishable when coming out the cavity. If we require the neutrons with 
relative velocity 5v z < 2 ■ Av z to be indistinguishable, the following relation must be satisfied 

< S «> 2< 2^' < 13) 
This gives a quantitative condition for two neutrons to be longitudinally indistinguishable 
throughout the cavity. At the same time, in the transection plane, the two neutrons would 
be projected into the same energy state ^0,0(^5 y) while approaching the exit of the cavity 
as the excited states are absorbed by the two upper absorbers, that is 

Equation (fl4l) is symmetric under the permutation of the two neutrons. Due to the Fermi- 
Dirac statistics of identical fermions, if two neutrons enter the cavity simultaneously and 
equation (fl3l) is satisfied, the spin wave function of the neutron pair must be antisymmetric, 
i.e. 

ri AB) = ^(1 + -)-!- +)) (15) 

when leaving the cavity. Inputting the values of h, m, and energy E 0t0 , E 0t i into equation 
( 1131) and choose T to be two orders greater than At (i.e. T = 10 2 • At), we get the critical 
value for Av z as 

A^ z < 0.6 ■ KT 3 ms" 1 . (16) 

Note that while choosing T = 10 ■ At, we then get Av z < 2 ■ 10 -3 ms _1 . Thus the neu- 
tron sources with a FWHM of velocity less than 10 -3 ms" 1 are capable of generating spin 
entanglement within our scheme. 
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Figure 3: A schematic setup of the devices for the generation of gravity induced entangled neutron 
pairs and the measurement of Bell inequalities. Here the arrows represent the motion of neutron, 
(a) is the energy filter of neutrons, including: two superconducting solenoid polarizers (SCS) and 
a spatial spin resonance (SSR) unit; (b) is the top view of the rectangular cavity (RC) of Figure [TJ 
(c) is a compound of two SCS acting as two separated neutron polarization analyzers. 



In the following we will present a practical experiment setup for generating entangled 
neutron pairs via the Earth's gravitational field using ultra cold neutrons (UCN, velocity 
v < 10 ms^ 1 [23] and we choose v = 5 ms" 1 for numerical evaluation hereafter). The three 
ingredients in our scheme which have been stated at the start of the paper include: the 
energy filter, the rectangular entangling cavity and the spin analyzers, shown in Figure |3j 



First, the monochromatic UCN of AA/A = Av/v < 1CT 3 (using values of equation fjl6H ) 
could be achieved via a Drabkin energy filter [24J which is composed of a polarizer, spatial 
spin resonance units, and a polarization analyzer. Using the superconducting solenoid- 
polarizer [22] (the left SCS in Figure. E(a)), the UCNs with different polarizations (parallel or 
antiparallel to magnetic field) have contrary behaviors: one polarization passes unhindered 
through the solenoid while the other one being reflected back. After this, the degree of 
polarization of UCN here can reach the level of 100% [26]. The spatial spin resonance units 
(SSR in Fig. E^a)) only flip the spin of neutrons for specific wavelength (velocity) [23]. The 
spin analyzer (the right SCS in Fig. Et^a)) selects the spin flipped neutrons which now have 
the monochromatic wavelength. A precision of AA/A ~ 10~ 2 has already been reached in 
for neutrons with A = 5 A. As the resolution of the A is inversely proportional to 
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the number of SSR units that neutrons pass through, an improvement to less than 10 3 is 
technically straight forward |28j . 

After the monochromatization, two neutrons are led to the rectangular cavity as Fig. [1] 
with relative angle of 9 ~ 10~ 2 (see Fig. 13(b)) in the horizontal plane ensuring their velocity 
in the transection plane is about 10 -2 ms _1 . The rectangular cavity with a length of about 
(5ms -1 ) • T would be enough for two neutrons coming out the cavity to be in spin singlet 
state. The out going neutrons can be guided to sufficient long distance as the depolarization 
effect of UCN in collision with materials is quite low, ~ 10~ 5 per collision [29]. Localized 
polarization analyzers can then be applied to the two well separated neutrons (see Fig. [3(c)). 
In particular, when two analyzers are placed in parallel, if one neutron passes through the 
solenoid polarizer then we can determine with certainty that the other neutron would be 
bounced backward from its solenoid polarizer though there is no interaction to project it 
to specific polarization state. A measurement of the Bell inequalities can be acquired by 
actively choosing different directions of the two solenoid polarizers separately. 

In conclusion, it is demonstrated that through a particular type of rectangular cavity 
two neutrons immersed in the earth's gravitational field will entangle with each other. This 
enables us to test the quantum nonlocality involving the gravity which is the only fundamen- 
tal interaction of nature having not yet been touched in practical entanglement generation 
so far. Due to the high detection efficiency of massive particle and the manageable large 
spatial separation between two neutrons (mean life of neutron at rest is r ~ 885.7s [30]), the 
proposed scheme in this work provides a simple and practical way for the implementation of 
nonlocality test of quantum entanglement and statistics in gravitational field, while a more 
conclusive test of local hidden variables theory would also be expected. 
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Appendix 

Separation of the variables 

Considering the potential of equation (J3J) in the region of x > 0, and y > 0, equation (jlj) 
can be expressed as 

+ W ) * M + + * * M = ^ (17) 

Defining ty{x,y) = ip(x)ip{y), the separation of variables, and E = E x + the above 
equation can then be expressed as 

+ ' 

These two equations can be solved in a similar way as equation (j2J). The solutions are 

^ n (x) = J\f n M(x/l - E n /e ) , tp m (y) =J\f m Ai(y/l - E m /e ) , 

where Ai are Airy functions; n, m G {0, 1,2, ■ ■ ■ }; A/" n is the normalization constant; Zo> £ o 
are the characteristic length and energy defined as Iq = h 2 ^ 3 / (y/2m 2 g) 1 ^ 3 , Sq = \Jmg 2 h 2 / '4. 
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Figure 4: The zeros of Airy function: Ai(— a n ) = 0, n G {0, 1, 2, • • ■ }. 

The eigenvalues of energy can be obtained by imposing the boundary condition i[) n (0) = 
Af n Ai(—E n /e ) = 0. We can get E n = a n e , where a n is the nth zero of Airy function, see 
Figure HI 
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Bouncing in the cavity 

In the transection plane of the cavity, the ground state in the gravitational potential is 
^o,o( x ,y) — V'oC^OV'o (y)- Inputting the numerical values E Q , J\f into the wave function we 
have 

25 

tt 0) o(s, V) ^ — Ai(x/6.59 - 2.34) • Ai(y/6.59 - 2.34) , (19) 

ol 

E 0fi = 2.22 • 10" 12 eV , (20) 

where x, y are in unit of /zm = 10 _6 m. In the ground state, (x) = (y) ~ 10.29 fim, 
(x 2 ) = (y 2 ) ~ 126.89 (/im) 2 , then Ax = Ay = 4.59 /un. Due to the Heisenberg uncertainty 
relation AxAp x > h/2 and p x = mv x , we have Av x > 6.86 ■ 10~ 3 ms _1 . The maximum 
velocity in x axis is v xmax = ^/2E /m ~ 1.46 • 10~ 2 ms _1 . The average velocity v x in x axis 
satisfies 



J xmax 



2 



< A« x . (21) 



T ■ v x 

The number of bouncing times of the neutron in the cavity is n = x , where l x is the 

Zlry 



"X 



side length of the transection plane. The uncertainty of the bouncing times n arises in 
regard of the variance Av x , and can be expressed as An = ^ x > n, which tells An > 1. 
Consequently, it is not distinguishable whether the neutron bounces in the cavity with odd 
or even number of times. 
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